www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/ybjzf_zuoye2.m
function result=ybjzf_zuoye2(a1,b1,a2,b2,mm) %a1是外重积分的下限 %b1是外重积分的上限 %a2是内重积分的下限 %b2是内重积分的上限 %积分函数cos(x)+sin(y)+3 %mm 是随机实验次数 xrandnum = unifrnd(a2,b2,1,mm); yrandnum = unifrnd(a1,b1,1,mm); fsum=sum(cos(xrandnum)+sin(yrandnum)+3); result=fsum*(b1-a1)*(b2-a2)/mm return ybjzf_zuoye2(0,4,1,3,10000) result = 24.4723